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Abstract: In the framework of a relativistic covari- 
ant Bethe-Salpeter model for the quark-antiquark system we 
present a renewed determination of the Cabbibo-Kobayashi- 
Maskawa matrix element \V c b\- Complementing an earlier 
analysis applied to the whole decay spectrum for B — > D'ev 
we now also employ the "zero-recoil method" that uses the 
end point of the decay spectrum (u = 1) and is suited for 
heavy-to-heavy transitions. The averaged experimental value 
extracted from the data at zero recoil, |V c i,|jT(o; = 1) = 
0.0343 ± 0.0015, then leads to \V cb \ = 0.0360 ± 0.0016. This 
value is somewhat larger than the one that uses the whole 
decay spectrum for the model analysis. We also contrast this 
result to a nonrelativistic model and to recent experiments on 
the B —* Dev semileptonic decay. 
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I. INTRODUCTION 

Within the standard model the extraction of the 
Cabbibo-Kobayashi-Maskawa (CKM) matrix elements 
\V c b\ (and \V u b\) is an outstanding topic of B-meson 
physics. Several ways have been utilized that are summa- 
rized, e.g., by the Particle Data Group Presently, the 
value of |T4h| extracted from inclusive decays is somewhat 
larger than from exclusive decays, e.g., in B — > D*ev. 

A fruitful method to extract \ V c b \ from exclusive decays 
is to reparameterize the decay data in such a way that 
they may be fitted by a smooth (e.g., linear) function 
|Vd,|.F(u;) of ui, where u = {m 2 B + m 2 D , — q 2 ) / {2mBra£,* ), 
and q 2 is the 4-momentum transfer. Doing so it is possi- 
ble to extrapolate to the point of zero recoil of the D*- 
meson, i.e., w = 1, that is not directly measurable. This 
procedure is particularly favored in the context of heavy 
quark expansion (HQET) || but also useful to compare 
to other approaches since in this context the notion of 
the whole decay spectrum is not needed to extract \V c b\- 
In HQET the value of T at zero recoil is normalized up to 
corrections of order (Aqcd/itt-c^) 2 (where m c ,b denotes 
the mass of the c- or &-quark). However the required fit- 
ting and extrapolation procedure leads to some errors, 
where the statistical error is under control and presently 
in the order of 5% jj) . 

Alternatively, quark models have been proven very use- 
ful as they provide not only predictions for T(u>) for all 
u) and \Vcb\, but also numerous testable results for quite 
different processes P-|l2[. A general overview on bound 
state models for heavy hadron decay form factors has 
been given by Ref . |l3| . Other approaches to the physics 
of heavy quark has profited from are QCD sum rules E3] 



and lattice QCD @. 

Since T{ui) is known for a quark model one may ask 
for the implications on the empirical value for \V c b\, if 
the zero recoil result is contrasted to the one obtained by 
using the whole decay spectrum. Both methods are fre- 
quently used but not yet compared to each other directly. 
In addition, relativistic quark models also allow us to de- 
scribe heavy-to-light transitions, in particular B — > 7r(p) 
important to determine \V u b\, see, e.g., Ref. Q. In this 
sense the heavy-to-heavy transitions provide an impor- 
tant test case and bear an impression of possible model 
uncertainties. 

In this context the merit of semileptonic B — > D^ev 
transitions may be considered twofold: As already men- 
tioned they provide an very good source to extract 
the Cabbibo-Kobayashi-Maskawa (CKM) matrix element 
\V c b\ that has to be contrasted to inclusive and nonlep- 
tonic decays. On the other hand, weak decays (in gen- 
eral) provide important complementary information for 
QCD-motivated modeling of the underlying quark struc- 
ture of mesons (in general hadrons). In addition, they 
may be considered useful to discuss the different rela- 
tivistic approximations used in this context. 

We choose an approach utilizing the instantaneous 
Bethe-Salpeter equation to treat the oq-system within 
a relativistically covariant formalism [p.7| . The model is 
able to describe the meson mass spectrum for low radial 
excitations. It has been applied to the calculation of lep- 
tonic decays, viz. decay constants, 77-decays |l8| , and to 
elastic form factors of mesons jlj| as well as to charmo- 
nium and bottomonium [ p0[ . Relativistic quark models 
have been investigated, e.g., in Refs. [pl|-p5|. 



II. THE BETHE-SALPETER APPROACH 

A. Solving the Bethe-Salpeter equation 

The Bethe-Salpeter approach provides a consistent 
treatment of two-body bound states as well as the cou- 
pling of an external field via the Mandelstam formal- 
ism p6| , p7| . In order to actually solve the bound state 
problem several reasonable approximations are necessary 
or practical: i) The quark propagators are assumed to 
be free propagators irrespective of confinement that is 
introduced via a confining kernel, ii) quark masses are 
assumed to be constant (i.e. constituent quark mass) 
which is reasonable for heavy quarks, since current quark 
masses and constituent quark masses needed for repro- 
ducing the mesonic mass spectrum are rather close to 
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each other, iii) we utilize ladder approximation for the 
interaction kernel, and iv) using an instantaneous inter- 
action in addition leads to computational advantages, as 
it provides RPA-type equations |28[ that can be solved 



f dp 
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of the Bethe-Salpeter equation, and it allows to substi- 
tute q 2 by -q 2 . The OGE kernel then reads 



by introducing an effective Hamiltonian 21 1 in a formally 
covariant way. The specific model used here, has been 
solved for the gg-systcm in|17]|l8| and applied to a wide 
range of phenomena [[19[|29| including the heavy quark 
sector ||[20) . Details of the model may therefore be found 
in the references given in the introduction. Here, I give 
a short survey and summarize some results. 

Within the approximations given above the po- 
integration in the Bethe-Salpeter (BS) equation may be 
performed. The resulting Salpeter amplitude in the rest 
frame of the bound state with mass M is the given by 



[^ G (p,p')<I>(p')] = V G ((p-p') 2 ) 



7°<fr( P ') 7° - 5 (7*(P ') 7 + (7*)*(P ') (7*) ) 



with the operator x = x/|x|, and 



V G (q 2 ) =tt 



4 a s (q 2 ) 
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where xp(p°,p) is the full Bethe-Salpeter amplitude. 
Note, that the relative momentum p = (joq, p) appearing 
in Eq. (|J) may be written in a covariant fashion [ pT|JT^ ] . 
The resulting Salpeter equation is then given by 



Here uji 

tion operators A^(p) = (a;, ± Zfj(p))/(2u;j) in obvious 
notation, where Hi(p) = 7°(7 • p + m,) is the standard 
Dirac Hamiltonian (for details see, e.g., Refs. [ pTp^ ] ). 

The dynamical input of the model is defined by a con- 
finement plus one gluon exchange (OGE) kernel, V = 
Vc + Vq. Confinement is introduced as a mixture of a 
scalar and a vector type kernel in the following way, 

[V c (p,p') $(p')] = Vg(( P - p') 2 ) [«(p') ~ 7° *(p')7°]. 

(3) 

Due to the instantaneous approximation it is possible to 
introduce the same spatial dependence (in the rest system 
of the meson) as used in the nonrelativistic case, viz. in 
co-ordinate space, 



where a s (q 2 ) is introduced as "running" coupling as dis- 
cussed in Refs. 

To solve the Salpeter equation numerically, Eq. (||) 
is rewritten as an eigenvalue problem (RPA-equations), 
see, e.g., Rcf. 021). This way it is possible to utilize the 
variational principle to find the respective bound states. 
To this end the Salpeter amplitude $ is expanded into 
a reasonable large number of basis states used as a test 
function. As a suitable choice of basis states we have 
taken Laguerre polynomials and found that about ten ba- 
sis states lead to sufficient accuracy, see also Refs. p7| , p^[ . 

For completeness, the parameters of the model given in 
Rcf. H are shown in Table [| These are the quark masses, 
the offset a c and slope b c of the confinement interaction 
Eq. (Q) and the saturation value a sa t = «s(q 2 0). 
( 2 )They are determined to give a good overall description 
of the meson mass spectrum (heavy and light mesons as 
well as charmonium and bottomonium) |P,|l7|- pO| , p9| . 



B. Current matrix elements 

Semileptonic decays are treated in current-current ap- 
proximation. For a transition b — > c the Lagrangian is 
given by 



Ccb = V cb K b j^ 



(7) 



with the CKM matrix element V c b and the Fermi con- 
stant Gf- The leptonic and hadronic currents Wti are 
defined by 



The mixture of a scalar and a vector spin structure has 
been introduced in order to give an improved description 
of the spin orbit splitting. Other mixtures have been ad- 
vocated in the literature, and also anomalous tensor-type 
confinement has been discussed, see, e.g., Ref. 0. How- 
ever, the consequences concerning, e.g., the mass spec- 
trum or Regge behavior have not been studied yet. 

For the OGE kernel, we chose the Coloumb gauge for 
the gluon propagator. This way it is possible to retain a 
covariant formulation within an instantaneous treatment 



3ti = huO- - 75^, 
The relevant transition amplitudes 



£>(*) | g y(l- 75 )&| B \j 



(8) 
(9) 



(10) 



for B — > D and B — > D* of the hadronic current can be 
decomposed due to the Lorentz structure of the current, 
thus introducing form factors. A standard representation 
of the form factors is given in terms of Fo(q 2 ), Fi(q 2 ) 
(or f+(q 2 ), f-(q 2 )) for 0~ — > 0~ transitions, and V(q 2 ), 
A (q 2 ), A 1 (q 2 ), A 2 (q 2 ) for 0" -> 1" transitions The 
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exact definitions, and further references have been given, 
e.g., in Ref. Q. Note that m| < q 2 < q 2 nax = (m B - 
m D (,)) 2 due to kinematical reasons. Helicity amplitudes 
H± and Hq in terms of the above form factors have been 
given by Korner and Schuler in a series of papers |30| and 
are compiled by the particle data group Q. Using the 
helicity amplitudes the formulas for the decay spectrum 
used here is given in Ref. (3^] and will no be repeated 
here. The respective decay rates into specific helicity 
states r±, Tq are also given in the literature, see, e.g., 
Ref. 0. 

To determine the form factors from the model, we fol- 
low the general prescription by Mandelstam j27j] , see, e.g. , 
Ref. |3j| for a textbook treatment. The lowest order 
contribution (relativistic impulse approximation) to the 
current (sometimes referred to as triangle graph) is given 
in Fig. |l|, and written as (consider, e.g., the anti-quark 
current, flavor indices suppressed) 



(D*,P D , \hM\B,P B ] 



d 4 p 



(11) 



MA", i---/ j (27r)4 

tr {fp D , (p - q/2) S$(P B /2+p- q) 
xV(l-76) 

xSf (Pb/2+p) T Pb ( P ) S^(-P B /2+p)} 

where p and p' denote the relative momenta of the incom- 
ing and outgoing qq pair, q = Pjj* —P B is the momentum 
transfer. The quark Feynman propagator is denoted by 
Sg . The Dirac coupling to point-like particles is consis- 
tent with the use of free quark propagators. In Eq. Jll| ) 
the amputated Bethe-Salpeter amplitude or vertex func- 
tion Tp(p) is given by 



(12) 



It may be computed in the rest frame from the equal time 
amplitude ^(p) using the Bethe-Salpeter equation 



= r( P ) = -i 



d 3 P ' 

(2tt)4 



[v( p ,p'Mp')] 



Xp(p)= S Ap X(Af,o)(A P 1 p) S A l- 



expansion Q. Following Ref. flj we for one introduce 
the ratios -Ri(w) and R 2 (oj), 

r 2 

1 



R 2 (u) = 



[m B + m D *) 2 



1 - 



<t 



(m B + m D ,) 2 



V(q 2 ) 
Ai(g 2 )' 

A 2 (q 2 ) 



(15) 
(16) 



where u> — [m B +m D » — g 2 )/(2mgmo«). 

For B — > D*ev decays the standard form factors may 
then be related to the ones used for the heavy quark 
expansion by J3lj 

„2 



Al(q 2 ) 

A 2 (q 2 ) 
V(q 2 ) 



K BD * 



1 



(m B + m D *) 2 

K B D* #2(")/l,4i 

k B d* Rii^hAxiu)- 



h Al (Lu), (17) 

(18) 
(19) 



where k B d' — (m B + mu»)/ (2y/m B rriD* )■ In the heavy 
quark mass limit (m c ,b ~ * °°) 

^xH-f^H, (20) 

iMw) 1. (21) 

where is a universal function known as Isgur-Wise 
function ||. 

For the case B — > Z?e^ the standard form factors are 



related to the heavy quark form factors via |31|, 



Finally, using the Lorentz transformation properties of 
the field operators that define the Bethe-Salpeter ampli- 
tude Q, we can calculate the Bethe-Salpeter amplitude 
in any reference frame via 



(14) 



where Ap is the pure Lorentz boost, and S Ap the corre- 
sponding transformation matrix for Dirac spinors. 

Due to the reconstruction of the full Bethe-Salpeter 
amplitude sketched above, the transition matrix element 
Eq. (O) is manifestly covariant. 



/+(<7 2 ) = ^ D h + (ui) - K BD h-(ui), (22) 

f-(<l 2 ) = K BD h -(u) - K BD h + (uj), (23) 

where = (m B ± mjj) / {2^Jni B mi)). In the heavy 

quark mass limit h+(uj) — > and h—(w) — > 0. 

In the model approach used here the heavy quark mass 
limit has been performed numerically by multiplying m c ^ 
with a large factor and keeping all other parameters as 
given in Table fl. To evaluate the transition matrix el- 
ements Eq. ([11 ) the meson amplitudes are then calcu- 
lated by diagonalizing the eigenvalue problem with the 
,_Jarge quark masses. Due to numerical reasons the heavy 
^ quark masses cannot be chosen too large. The function 
resulting from this numerical limiting procedure is then 
defined to be the Isgur Wise function £(w) of the Bethe- 
Salpeter-model, where £(1) = 1.00 within 0.1%. This 
function is shown as a solid line in Fig. ||. Note that at 
this stage £ does not include radiative corrections that 
will be given below. In the same fashion the ratios R± 
and R 2 tend to unity within less than 0.1% when numer- 
ically increasing the heavy quark masses. 

For finite masses the experimental ratios R\^ 2 assum- 
ing constant values have recently been extracted by 
CLEO 101. The latest values are I02L 



Rx = 1.24 ±0.26 ±0.12, 
R 2 = 0.72 ±0.18 ±0.07, 



(24) 
(25) 



C. Form factors 

The analysis of experimental data on heavy-to-heavy 
transitions now widely uses the notion of heavy quark 



that have to be contrasted to the long-dashed (Ri) and 
dashed-dotted (R 2 ) curve shown in Fig. g. The model 
ratios vary slowly by roughly 10% over the whole u> range, 
which is smaller than the experimental error. 
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III. RESULTS 

Utilizing the Bethe-Salpeter model to describe mesons 
as qq states the exclusive decay spectra for B — ► D*ev 
and B — > Dev have been calculated and compared to the 
experimental data. Earlier the CKM matrix element \ V c b\ 
has been determined by a least squared fit to the whole 
spectrum of B — > D*ev ||. We have now redone this 
analysis on the basis of the improved data and also in- 
cluded the recently measured B — ► Dev decay spectrum. 
In addition, we present a new analysis for this model ap- 
proach at zero recoil of B — > D*ev that is commonly 
used for heavy-to- heavy transitions to extract | V cb | . This 
enables us to compare the difference between the energy 
dependent and the zero-recoil analysis quantitatively, at 
least for the Bethe-Salpeter model discussed here and the 
nonrelativistic approach given earlier |33] . 



A. B — > D*eu decay 



We now turn to the extraction the CKM matrix ele- 



ment | V^b 
given by 



(II 



D* 



The differential decay rate for B — > D*ev is 



- m D ,) 2 

f\V cb \ 2 Tl 
% + m%. - 2um B m D . 



r 2 

m D ,(rn B 



48tt 3 

f 



4w 



1 



(m B - m D * 



(26) 



The formula is written in a way that Tv (w) reduces 
to Tr3*{w) — > r)A£(w) = m the heavy quark mass 
limit, where tja denotes the radiative corrections pl| , p6[ . 
For finite masses the function To* (w) contains all the 
symmetry breaking effects. 

Experiments are given in a way that all well known 
factors are divided out in the decay rate and only 
I V c b\ 3~d* (w) is left over. The corresponding data points 
of a recent CLEO measurement are shown in Fig. 0. The 
result is particularly smooth and may be fitted by a lin- 
ear curve. The fit to the data done by the CLEO col- 
laboration is also shown in Fig. || as dashed-dotted line. 
Other lines reflect the model results utilizing different as- 
sumptions. The solid line is calculated using the exact 
formula for the decay rate as given, e.g., in fl]|3Cj (i.e., 
with i?i ; 2 w-dependent) divided by the same factor as 
the experiments are. The CKM matrix element \V c b\ is 
then determined by a least squared fit to all data points. 
Radiative corrections have been included in the domi- 
nant form factor /i j 4 1 (oj), expressed as an overall factor 
rj a |5lj • Unlike earlier estimates that imply a correction 
of approx. 1% J3l| (i.e., smaller than the model uncer- 
tainty and therefore neglected in the earlier analysis 
a recent two loop calculation leads to substantial value 
of T] A = 0.960 ± 0.007 @, which has to be included into 
the analysis. The result is \V cb \ = 0.0339 ± 0.0010. To 



see the model dependence of the different analyses used 
in this context we now take Fd»(<*>) = rij^hj^iu)), R\ 
and i?2 constant, i.e., i?i(l), #2(1), and this CKM ma- 
trix element that leads to the long-dashed line shown in 
Fig. H| It is obvious that the curve slightly deviates from 
the solid one that includes the w-dependence of R\ and 
i?2- For the same form factor however using the value 
of \V c b\TD*{u = 1) from the CLEO fit is given by the 
short-dashed curve. The resulting CKM matrix element 
is obviously larger by ~ 8%. 

The function Td* (uj) that leads to Fig. |^ can be ap- 
proximated by a quadratic fit to 

^}(cj) = Tn* (1) (1 - P \ (u-l) + c{oj - l) 2 ) , (27) 

with the parameters p A and c. The slope of Tr>* (uj) 
extracted by CLEO g assuming a linear dependence on 
oj (c = 0) is p\ = 0.84±0.13±0.08, and show as dashed- 
dotted line in Fig. |^. The respective parameters for the 
quadratic fit of the curve discussed are shown in Tab. ||. 

Within the notion of the heavy quark effective theory 
the form factor — 1) can be expanded into or- 

ders of ^QCD/n^c,b- Since symmetry breaking effects in 
semileptonic decays are of second order only J37[ , lowest 
order terms are usually written as 



T d *{uj = 1) = tja{1 +$i/ m i), 



(28) 



where Si/ m 2 has to be determined. Corrections vary from 
-S 1/m 2 = (3 ±2)% §1|§ to -5 1/m * = (5.5 ±2.5)% g§. 
A recent discussion and appreciation of the different ap- 
proaches is given by Martinelli j4fj]. The relativistic 
model discussed here leads to value of —6i/ m 2 — 0.5%. 

The values for Td'{u = 1)|V*| that have been ex- 
tracted by different experiments are quite consistent and 
lead to an overall fit of |4l| , j32"f 



T d ,(oj = l)\V cb \ = 0.0343 ±0.0015 



(29) 



for the empirical slope parameter given in the last line 
of Tab. ||. Using this value and the radiative corrections 
given above we extract the CKM matrix element \V cb \ for 
the relativistic Bethe Salpeter model to be 

\V cb \ = 0.0360 ±0.0016 zero-recoil, (30) 
\V cb \ = 0.0339 ± 0.0010 full spectrum. (31) 

This is the main result that shows the potential model 
dependence of the zero-recoil method that adds to the 
statistical uncertainty. A similar renewed analysis for 
the nonrelativistic model given before now leads to 
\Vcb\ = 0.037 ± 0.002 that is larger by 7% compared to 
the spectrum dependent analysis. 



B. B — » Dev decay 
The differential decay rate for the B — * Dev is given 

by & 
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x(w-l) 3/a |Vc&| 3 *£H (32) 

where again Td{oj) reduces to the Isgur Wise function 
in the heavy quark mass limit. Recent experimental 
data for the relevant part of the spectrum | V c b\ Fd{w) are 
shown in Fig. [|. The linear fit to the data given by the 
CLEO collaboration is also shown (as a dashed-dotted 
line). The model results utilizing the full spectrum de- 
pendent analysis leads to the solid line in Fig. [|. For 
comparison the result of the zero recoil method utilized 
in the previous paragraph is also shown. The radiative 
corrections have been assumed to be in the same order as 
in the B — > D*ev transition. Obviously the model is ca- 
pable to provide a good description of the experimental 
data. 



the statistical error. This result may also be relevant for 
other more "model independent" analyses. 

For comparison the B — > Dev decay has been calcu- 
lated for the different assumptions and also leads to a 
good overall description of the experimental spectrum. 
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IV. SUMMARY AND CONCLUSION 

We have analyzed the exclusive decay rates of B — > 
D*ev and B — > Dev within a relativistic constituent 
quark model. The interaction kernel has been taken in- 
stantaneous. This way the Bethe Salpeter equation re- 
duces to a Salpeter equation as given in Eq. (||). The 
interaction consists of a one gluon exchange evaluated in 
the Coulomb gauge and a linear confinement given in co- 
ordinate space. The model parameters have been fixed to 
describe the mass spectrum of all observed mesons (not 
only heavy mesons) in a satisfactory manner Jl7|-|2C|l . The 
interaction current to describe the weak decay process 
has been introduced via the Mandelstam formalism. To 
this end the (instantaneous) amplitude has been recon- 
structed using the Lorentz transformation properties of 
the field operators. 

The only parameter left to describe the exclusive decay 
spectra is the CKM matrix element \V c b\. For B — ► D*ev 
two methods have been compared. One uses the complete 
spectrum, viz. the functional dependence of Try* , emerg- 
ing from the quark model. The CKM matrix element 
is then fixed by a least squared fit. The other analysis 
utilizes the "zero recoil" method used in the context of 
heavy quark expansion. Here only the empirical value of 
IKhl^D* (1) is used that is gained from an extrapolation 
of the experimental data (e.g. by a linear fit) to the zero 
recoil point u> = 1 and R2 assumed constant. The 
CKM matrix element \V c b\ is then extracted using a sin- 
gular value of the model at Tv (1)- Clearly, this method 
is not consistent with the underlying model, but widely 
used to extract \V c b\ from the zero recoil point. Com- 
paring the full model result to the zero-recoil result leads 
to different values for the CKM matrix elements \V c b\ by 
approximately 8%. In view of this result it seems obvi- 
ous that this kind of model uncertainty steaming from 
the different treatment of the w-dependence of the spec- 
trum may show up in the determination of \V c b\ besides 
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TABLE I. Parameters of the BS model. 



m Ut d m s m c m b a c b c r a sa t 

[GeV] [GeV] [GeV] [GeV] [GeV] [GeV/fm] [fm] 

0.200 0.440 1.738 5.110 -1.027 L700 0.1 0.391 



TABLE II. Parameters for a quadratic fit to the CLEO 
data and the BS model. First line in heavy quark mass limit, 
other lines for physical masses. 



p\ 


c 


model 


0.83 


0.34 


for £bs(w) 


0.76 


0.30 


for h Al (w) 


0.69 


0.34 


for T{lu) 


0.92 ± 0.64 ± 0.40 


0.15 ± 1.24 ±0.90 
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FIG. 1. Pictorial demonstration of one particle approx- 
imation of the irreducible interaction kernel, filled circle de- 
notes the inverse quark propagator 
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FIG. 2. Form factors c$i (long-dashed), R2 

(dashed-dotted), and hA ± for the B — > D* transition. The 
Isgur-Wise function £ of heavy quark mass limit is shown as 
solid line. Radiativ corrections not included here. 
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FIG. 3. V c i,|JT(tj) as a function oil xj for B — > D*ev decays. 
Data points from CLEO and their linear fit (with c = 0, 
dashed-dotted line). Model best fit is displayed as solid line. 
Long-dashed with the same value for \V c b\ at u = 1 but Ri t 2 
constant. Dashed line uses V c i from CLEO analysis and Ri t 2 
constant (i.e. zero-recoil method) for comparison. 
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FIG. 4. I V c b |^"(w) as a function oY lj for B Dev decays. 
Data points from CLEO and their linear fit (with c = 0, 
dashed-dotted line). Model best fit is displayed as solid line, 
same V c b as in Fig. [| 
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